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Abstract 

The  objective  of  this  program  is  to  investigate  approaches  which  make  use  of  several 
classes  of  electroactive  materials  developed  recently  for  harvesting  electric  energy  during 
walking.  These  newly  developed  electroactive  materials  including  the  electrostrictive 
PVDF  based  polymers  and  piezoelectric  relaxor  single  crystals  have  shown  order  of 
magnitude  improvement  in  terms  of  the  elastic  energy  density  in  comparison  with  the 
traditional  piezoelectric  materials  and  much  improved  the  electromechanical  coupling 
factors,  both  of  which  are  required  for  high  performance  energy  harvesting  systems. 
Furthermore,  in  this  program,  through  a  systematically  analysis,  we  have  shown  that  the 
input  mechanical  energy  density  to  an  electroactive  material  and  the  electromechanical 
conversion  efficiency  of  a  properly  designed  energy  harvesting  system,  i.e.,  by  properly 
designed  electric  controlling  circuits,  can  be  much  higher  than  these  of  the  electroactive 
material  itself.  Using  this  principle,  we  demonstrated  that  an  electric  energy  output  of  39 
mJ/cc  with  a  10%  efficiency  can  be  obtained  from  the  electrostrictive  PVDF  polymers. 
Using  1-3  composites  with  ferroelectric  relaxor  single  crystals,  harvested  electric  power 
density  of  96  mW/cc  has  been  obtained  (at  a  frequency  of  4  Hz).  The  experimental 
results  have  revealed  several  issues  in  utilizing  these  electroactive  materials  for  the 
electric  energy  harvesting:  (1)  The  dielectric  hysteresis  loss  at  high  electric  fields  for  the 
electrostrictive  PVDF  polymers  should  be  avoided  in  order  to  achieve  high  efficiency;  (2) 
For  the  1-3  composites  with  the  ferroelectric  relaxor  single  crystals,  a  DC  bias  field  is 
needed  in  order  to  stabilize  the  material  in  the  piezoelectric  state;  (3)  Because  of  the 
limitation  of  the  commercially  available  power  electronic  components  as  well  as 
relatively  small  electric  energy  involved  in  the  energy  harvesting  systems  (~  1  watts  or 
below),  it  is  preferred  to  operate  the  energy  harvesting  system  to  below  1000  volts  in 
using  active  control  of  the  energy  harvesting  system  with  the  commercially  available 
power  electronic  components. 


I.  Introduction 

The  modem  warfare  is  marked  by  sophisticated  electronics  used  in  almost  all  parts  of  the 
weaponry  systems  which  enhance  significantly  the  military  capability  and  reduce  the 
battlefield  causality.  One  of  the  prices  for  this  increased  reliance  on  the  electronics  is  the 
increased  weight  of  the  batteries.  Carrying  these  batteries  can  increase  the  pack  weight  to 
upwards  of  80  lbs.  Our  goal  in  this  study  is  to  develop  a  device  that  “harvests”  some  of 
the  mechanical  energy  during  walking  and  converts  it  to  electrical  energy  to  charge  the 
batteries  so  that  the  amount  of  batteries  that  must  be  carried  for  each  day  in  the  field  can 
be  reduced  significantly.  Furthermore,  through  this  investigation,  we  also  intend  to 
establish  general  guidelines  for  the  energy  harvesting  using  the  electroactive  materials  to 
achieve  the  highest  possible  efficiency  and  maximum  possible  electric  power  output. 

In  a  typical  energy  harvesting  process,  an  electroactive  material  absorbs  the  external 
mechanical  energy,  and  through  a  proper  electromechanical  conversion  process  to 
convert  that  energy  into  electric  form.  Finally,  this  converted  electric  power,  through  a 
power  electronic  interface,  is  delivered  to  the  load.  Therefore,  several  steps  must  be 
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taken  in  designing  an  energy  harvesting  system  to  effectively  convert  the  mechanical 
power  available  in  an  external  medium  into  the  electric  power  used  by  an  electric  load: 

(i)  mechanical  impedance  match  of  the  electroactive  material  system  to  the  environment 
to  realize  maximum  mechanical  energy  transfer  into  the  electroactive  material  system; 

(ii)  high  electromechanical  conversion  efficiency  of  the  electroactive  material  system; 

(iii)  a  proper  power  electronics  to  transfer  most  effectively  converted  electric  power  into 
the  electric  load.  Furthermore,  we  use  the  electroactive  material  system  which  couples 
the  electroactive  material  with  a  properly  designed  electronics  to  “smartly”  tune  the 
material  parameters  during  the  energy  harvesting  cycle  for  the  energy  harvesting  system. 
These  electronic  additions  can  improve  the  electromechanical  conversion  efficiency 
and/or  the  mechanical  impedance  match  to  the  external  mechanical  environment.  That  is, 
through  properly  designed  power  electronic  circuit,  the  energy  conversion  efficiency  of 
an  electroactive  material  system  can  be  much  higher  than  that  of  the  electroactive 
material. 

II.  Results 

(i)  Optimizing  the  energy  harvesting  cycle  to  achieve  higher  electric  energy  output 
and  efficiency/11  Although  the  electromechanical  conversion  efficiency  for  a  given 
electroactive  material  is  a  fixed  value,  for  energy  harvesting  from  a  given  mechanical 
environment,  the  amount  of  mechanical  energy  as  well  as  the  system  efficiency  can  be 
improved  markedly  by  optimizing  the  electric  boundary  conditions  during  the  energy 
harvesting  cycle.  For  example,  presented  in  Table  I  is  a  comparison  of  the  maximum 
electric  energy  density  harvested  Wimax  and  system  effective  coupling  factor  k  for  a 
given  electroactive  polymer  (electrostrictive  PVDF)  under  a  maximum  150  MV/m  field 
and  20  MPa  stress.  Apparently,  the  electric  power  output  depends  crucial  on  the 
operation  conditions  of  the  electronic  controlling  circuits  and  can  vary  over  a  large 
range. 

Table  I. 


Electrical  Boundary  Conditions 

W1Max(J/cc) 

k 

Cycle  1  -  Constant  field  and  open-circuit 

0.14 

0.523 

Cycle  2  -  Constant  field 

0.99 

0.815 

Cycle  3  -  Open  circuit 

0.826 

0.834 

Cycle  4  -  Passive  diode  circuit 

0.082 

0.392 

Taking  electroactive  polymers  as  an  example.  For  mechanical-to-electrical  energy 
harvesting,  one  can  apply  various  electrical  boundary  conditions  during  the  stress 
variation  cycle  to  break  the  symmetry  of  stretching  and  contracting  of  the 
electrostrictive  polymer,  and  thereby  harvest  energy.  In  order  to  maximize  the 
harvesting  energy  density,  the  boundary  conditions  have  to  be  designed  so  that  the  area 
enclosed  in  the  loop  (see  figure  1)  is  as  large  as  possible  without  exceeding  the 
limitations  of  the  polymer  material  (e.g.,  maximum  stress,  breakdown  field). 
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Figure  1.  Energy  harvesting  cycle  where  the  electric 
energy  harvested  is  the  area  enclosed  by  the  strain- 
stress  (S-T)  or  charge  density-field  (D-E)  loop  (the  two 
are  the  same) 


Ideally  the  energy  harvesting  cycle  should  consist  of  the  largest  loop  possible,  bounded 
only  by  the  limitations  of  the  material.  However,  actual  implementation  of  the  optimal 
energy  harvesting  cycle  in  an  energy  harvesting  device  may  be  difficult  to  achieve.  The 
4  different  energy  harvesting  cycles  listed  in  Table  I  correspond  to  the  electrical 
boundary  conditions  that  can  be  applied  to  the  device  fairly  easily  with  power 
electronic  circuitry.  The  following  is  a  summary  of  the  analysis  for  these  four  cycles 
using  an  electrostrictive  material  under  a  DC  bias  field.  The  results  can  be  easily 
extended  to  the  piezoelectric  material. 

Electrostriction  is  generally  defined  as  quadratic  coupling  between  strain  and  electrical 
field.  The  strain  Sy  and  the  electric  flux  density  Dm  are  expressed  as  independent 

variables  of  the  electric  field  intensity  En  and  the  stress  Ty  by  the  constitutive  relation  as: 


S..  =  s  T..+M  ..E  E  , 

tj  ijkl  y  mmj  n  m’ 

D  =e  E  +2M  ..E  T.., 

m  mn  n  mmj  rt  ij 


(i) 


Where  sfjkl  is  the  elastic  compliance,  Mnmij  is  known  as  the  electric-field-related 

electrostriction  coefficient,  and  emnis  the  linear  dielectric  permittivity.  An  isotropic 

electrostrictive  polymer  film  contracts  along  the  thickness  direction  and  expands  along 
the  film  direction  when  an  electric  field  is  applied  across  the  thickness.  Assuming  the 
only  non-zero  stress  is  applied  along  x, ,  the  constitutive  relation  is  then  simplified  as: 


where 


s  =  sT  +  me2 
D  =  sE  +  2MET 


(2) 
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E  =  E3,D  =  D3,S  =  sv 
T  =  Tvs  =  sn,M  =M3v 
e  =  £33 

The  electromechanical  coupling  factor  is  defined  as: 

k= 

yw;+w2’ 

where  Wx  +  W2  is  the  input  mechanical  energy  density,  and  Wl  is  the  output  electrical 

energy  density.  The  coupling  factor  is  often  associated  with  a  specific  set  of  electrical 
boundary  conditions,  as  will  be  discussed  in  the  following.  In  all  the  discussions,  the 
application  of  stress  cycle  starts  from  zero  to  a  maximum  value  TMax,  and  then  returns 
to  zero. 


•  Energy  harvesting  cycle  #1  -  Constant  field  and  open-circuit  electrical  boundary 
conditions 


Figure  2.  Energy  Harvesting  Cycle  under  Constant  Field  (l- 
2)  and  Open  Circuit  (2-3)  Electrical  Boundary  Conditions 


In  this  cycle,  The  material  is  stress  free  at  state  1  of  the  cycle,  and  an  electric  field  £0is 
applied  and  kept  constant  as  the  stress  is  increased  to  Tmm ,  ending  in  state  2.  The 

electrostrictive  device  is  then  open-circuited  when  the  stress  is  removed,  ending  in  state 
3.  From  state  3  the  electric  field  E0  is  re-established  in  the  material,  returning  to  the 
original  state  1.  The  total  mechanical  input  energy  density  available  for  the  energy 
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harvesting  for  such  a  cycle  can  be  shown  to  be  Wx  +W2  =  —  sT^ 


2 

max 


while  the  energy  W, 


can  be  calculated  as  follows. 


As  the  material  is  in  open-circuit  condition  from  state  2  to  3,  the  electric  flux  density  is 
constant  and  is  given  by 
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D  —  £Eq  +  2ME0Tmax 


(3) 


The  electric  field-intensity  as  the  stress  is  re'moved  is  therefore  given  by: 

E  =  E  - — 

0  e+2  MT 

This  field-intensity  reaches  its  peak  value  E  when  the  applied  stress  becomes  zero,  or 


(4) 


f  2  MT 

EP=E0  1+ - =-  =E,{l+y),  (4) 

V  £  J 


Where  y  corresponds  to  the  relative  change  in  dielectric  constant  due  to  applied  stress, 
and  is  defined  as  follows. 


2MT 

y  ~~ _ max 


(5) 


For  existing  electrostrictive  materials,  y  is  less  than  1.  Now: 

W2  =  -jVdS  -2%-M'T^ 

=\{xL-eyX) 


(6) 


(7) 


The  coupling  factor  is  therefore  given  by 


k  =  IMEq 


-\fs£ 


(8) 


The  electromechanical  coupling  factor  k  for  this  cycle  is  that  used  in  the  literature  as 
the  coupling  factor  for  the  material  (except  that  Eo  in  here  is  the  lowest  field  during  the 
cycle).  If  E0  is  chosen  so  that  the  peak  electric  field  intensity  Ep  is  the  maximum 

allowable  due  to  material  constraints  Emax ,  the  maximum  harvesting  energy  density  is 
given  by 


W  = 

^lmax 


r 


(l+y) 


■MTE; 


max — max 


(9) 


And  the  coupling  factor  associated  with  this  energy  density  is  given  by 
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k  - 


(10) 


•  Energy  harvesting  cycle  #2  -  Constant-field  boundary  conditions  during 
stressing  and  unstressing  of  material 


Another  possible  method  of  imposing  electrical  boundary  conditions  is  to  keep  the 
electric  field  constant  as  the  material  is  stressed,  change  the  field  to  a  different  constant 
value,  and  then  remove  the  stress.  Such  a  cycle  is  shown  in  Figure  3. 

A  constant  electric  field  E0  is  applied  from  state  1  to  state  2  as  the  stress  is  increased 

toT^.  From  state  2  to  state  3  the  electric  field  is  increased  from  £0to  E1 ,  and  then 
kept  constant  as  till  the  stress  is  reduce  from  7^  to  0  from  state  3  to  state  4.  Finally,  at 
zero  stress  the  electric  field  is  reduced  to£0,  returning  to  state  1.  In  the  mechanical 

stress/strain  plot  the  path  1-4  and  2-3  are  not  parallel,  which  is  due  to  the  stress 
dependence  of  dielectric  constant.  The  converted  energy  can  be  shown  to  be 


W.  =  T  M 
1  max 

1 


(£?-£o2) 


(11) 


The  input  energy  density  W2  =~sTfa ,  and  the  coupling  factor  is  therefore  given  by 


k  = 


i 


(12) 


The  maximum  energy  harvesting  density  and  coupling  occurs  when  E0  is  set  to  be  zero 
and  El  is  set  to  Emm . 


Figure  3.  Energy  Harvesting  Cycle  Under  Constant 
Electrical  Field  Conditions  as  the  Material  is  Stressed 
and  Unstressed 

•  Energy  harvesting  cycle  #3  -  Open-circuit  boundary  conditions  during  stressing 
and  unstressing  of  material 

Another  type  of  electrical  boundary  condition  is  to  have  open-circuit  boundary 
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conditions  as  the  stress  is  applied  and  removed,  as  shown  in  Figure  4.  From  state  1  to 
state  2  and  state  3  to  state  4,  the  electric  field  changes  automatically  as  the  stress 
changes.  From  state  2  to  state  3  and  from  state  4  to  state  1,  the  electric  field  is  changed 
through  the  electrical  interface.  We  define  the  field  at  state  1  to  be  E0  and  the  field  at 
state  3  to  be  E4 . 


And  the  energy  harvesting  density  is 


W1  = 


MT 


1  +  y 


■[El -El] 


(13) 


Its  associated  coupling  factor  is  given  by 
k  = 


(1  +  y)M[E2-El) 


(l  +  yfsTmax+(l  +  y)M[E2  - El]-2M2Tmsy.E2  / £ 


2  v  •'  niiiA  v-  / /- -  l— '4  "u  j  A  max  4 

The  maximum  harvesting  density  then  occurs  then  E0  =  0  and£4  =  Emax . 


(14) 


Figure  4.  Energy  Harvesting  Cycle  Under  Open-Circuit 
Conditions  as  Stress  is  Applied  and  Removed 

•  Energy  Harvesting  Cycle  #4  -  Passive  diode  circuit  for  energy  harvesting 

A  circuit  that  has  been  proposed  for  electrostatic-based  energy  harvesting  can  also  be 
used  with  electrostrictive  materials,  and  is  shown  in  Figure  5.  The  circuit  uses  high 
voltage  diodes  for  passive  switching.  The  main  advantage  of  this  circuit  design  is 
simplicity  (This  is  the  circuit  used  in  the  SRI  International  energy  harvesting 
experiment). [2] 


8 


Figure  5.  Passive  Diode  Circuit  for  Energy 
Harvesting 


The  energy  harvesting  cycle  is  shown  in  6.  However,  as  the  circuit  is  passive,  the 
voltage  change  across  the  device  occurs  only  due  to  the  electrostrictive  effect.  As  a 
result,  it  can  be  shown  that  the  voltages  VL  and  VH  are  constrained  by  the  following 
condition: 

vL<vH<yvL  as) 

These  constraints  severely  impair  the  energy  harvesting  density  and  coupling  factor 
when  y  is  small. 


T  E 

Figure  6.  Energy  Harvesting  Cycle  for  Passive  Diode 
Circuit 


(16) 


(17) 
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And  it  associated  coupling  factor  is 


k  = 


YMEj 


]j(i +r)sTma+ME2m 


r 

2(1 +r) 


(18) 


As  can  be  seen  from  Table  I,  the  passive  diode  circuit  does  not  provide  high 
energy  density  output  and  efficiency.  In  an  early  SRI  experiment,  it  was  claimed  that 
a  high  output  electric  energy  density  was  harvested.  The  result  is  not  correct  and  the 
error  is  due  to  the  fact  that  they  have  reported  the  discharge  electric  energy  to  Vh  and 
have  not  included  the  energy  input  from  Vi  to  charge  the  electroactive  polymer  which 
is  very  significant. 

(ii)  Experimental  results  of  electric  energy  harvested  using  the  electrostrctive  PVDF 
polymer.  As  shown  in  Table  I,  one  of  the  energy  harvesting  cycle  which  can  lead  to  a 
high  energy  density  and  efficiency  is  the  cycle  3.  In  this  energy  harvesting  cycle,  a 
constant  DC  electric  field  is  applied  when  the  material  is  stress  and  unstressed.  Because 
of  the  limitation  of  the  power  electronics,  the  operation  voltage  of  the  energy  harvesting 
system  is  preferred  to  be  below  1000  volts.  For  the  electrostrictive  PVDF  polymers 
which  require  an  operation  field  up  to  100  MV/m,  1000  volts  limit  will  require  the 
thickness  of  each  electrostrictive  polymer  film  to  at  10  pm  and  a  multilayer  approach 
needs  to  be  taken  for  the  device.  Because  of  the  earlier  termination  of  the  program,  we 
were  not  able  to  develop  a  power  electronic  circuit  as  well  as  the  multilayer  device  with 
each  layer  thickness  at  10  pm  and  to  perform  the  experiment.  To  quickly  demonstrate 
the  concept  of  utilizing  smart  electronics  to  tune  the  electroactive  material  parameters 
during  the  energy  harvesting  cycle  to  achieve  higher  efficiency  and  higher  input 
mechanical  energy  density,  we  approximated  such  an  electric  condition  with  a 
sinusoidal  electric  signal  (which  eliminates  the  limitation  of  the  power  electronics). 
Detailed  analysis  indicates  that  when  the  AC  stress  signal  and  AC  electric  signal  are  out 
of  phase  by  90°,  the  system  reaches  highest  energy  density  harvested.  The  experimental 
results  with  electrostrictive  PVDF  confirmed  this  prediction.  As  shown  in  Table  n,  an 
energy  density  of  40  mJ/cm3  can  be  achieved  with  the  electrostrictive  PVDF  polymer 
(corresponding  to  a  coupling  factor  k=0.31).  In  the  experiment,  the  maximum  electric 
field  is  67  MV/m  (much  lower  than  150  MV/m  used  in  the  calculation  of  Table  I,  which 
explains  why  a  lowered  energy  density  and  efficiency  obtained  here  compared  with  that 
predicted  in  Table  I).  This  lower  field  is  due  to  the  limitation  of  the  electric  breakdown 
of  the  test  specimen  as  well  as  the  increase  dielectric  loss  at  higher  fields.  Even  with 
these  limitations,  the  electric  energy  density  and  efficiency  are  already  much  higher 
than  all  the  previously  reported  values  (although  SRI  claimed  to  have  higher  energy 
density,1 [2]  their  value  is  not  that  from  the  whole  stress  loop  and  therefore,  is  not  the  true 
energy  density  harvested.  This  has  been  analyzed  in  detail  in  (i))[1]. 
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Table  II 


Material 

Specific  energy  (mJ/cc) 

PVDF* 

0.5 

0.044 

PZT* 

1.5-5 

2.1 

Electrostrictive 
PVDF  polymer 

10 

39.4 

*MIT,  IEEE  conference  on  wearable  computing  1998 


(iii)  Development  of  1-3  single  crystal-polymer  composites  with  high  electromechanical 
conversion  efficiency. 

In  our  proposal,  two  classes  of  electroactive  materials  were  proposed  for  use  for  the 
energy  harvesting:  the  electrostrictive  PVDF  polymer  and  1-3  piezocomposite  using  the 
high  efficiency  piezoelectric  single  crystals.  In  this  program,  the  fabrication  process  of  1- 
3  single  crystal  (PMN-PT)-polymer  composite  was  investigated  in  order  to  achieve  high 
coupling  factor.  Single  crystals  of  0.67PMN-0.33PT  oriented  along  <001>  direction  were 
used  in  the  investigation  because  of  the  high  coupling  factor  (>0.9).  The  single  crystals 
were  poled  at  room  temperature  under  a  field  of  10  kV/cm.  In  the  1-3  composite 
fabrication  process,  extra-precaution  has  to  be  taken  in  order  not  to  cause  damage  to  the 
single  crystals.  The  poled  single  crystals  of  1  mm  thick  (<001  >  oriented)  were  first  cut 
along  the  <100>  direction  using  the  Automatic  Dicing  Saw  with  the  a  curf  width  of  100 
pm.  The  cutting  process  was  followed  by  filling  of  Mereco  1650  series  epoxy  to  form  a 
2-2  composite.  The  2-2  composites  were  cut  in  perpendicular  direction  and  then  filled 
with  the  epoxy  to  form  the  final  1-3  composites  (see  figure  7).  The  1-3  composites  were 
polished  to  a  thickness  of  about  0.7  mm  (see  figure  7(a),  thickness  along  the  z-direction) 
and  sputtered  Au  films  of  20  nm  thick  were  used  as  the  electrodes.  The  composites  were 
poled  along  the  thickness  direction  in  silicon  oil  for  5  min  under  a  DC  field  of  1  kV/cm. 
Several  1-3  composites  with  piezoelectric  single  crystal  PMN-PT  were  fabricated  and 
characterized.  Shown  in  figure  7(b)  is  a  top  view  of  a  1-3  single  crystal-polymer 
composite  fabricated,  where  the  dark  areas  are  single  crystals  and  the  less  dark  areas  are 
epoxy. 


PMN-33PT  (phase  I)  Epoxy  (phase  2) 
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Impedance 


1x10*  2x10*  3x10*  4x10s 

Frequency(Hz) 


High  electromechanical  coupling 
factor  was  achieved  in  these  composites. 
Shown  in  figure  8  is  the  electric 
impedance  curve  for  a  1-3  single  crystal- 
polymer  composite.  The  thickness 
coupling  factor  kt  of  1-3  single  crystal- 
polymer  composite  deduced  from  the  data 
is  0.85.  Figure  9  shows  how  kt  varies  with 
the  volume  content  of  the  single  crystal  in 
the  1-3  composites.  In  the  volume  content 
from  0.4  to  0.6  range,  kt=0.85.  In 


Figure  8.  Electric  impedance  curve  for  a  1-3 
single  crystal-polymer  composite  where  the 
solid  curve  is  the  phase  angle  and  the  curve 
with  dots  is  the  absolute  value  of  electric 
impedance 

comparison,  1-3  composites  with 
piezoceramics  PZT  show  a  coupling 
factor  of  0.59.  This  difference  is 
significant  since  the  energy  conversion 
efficiency  is  proportional  to  the  square 
of  the  coupling  factor  (72%  for  single 
crystal  composites  vs.  35%  for  PZT 
composites).  Moreover,  as  has  been 
shown  in  the  results  from  (i),  this 
difference  in  the  coupling  factor  can  also 


0.4  0.6  0.8  1.0 

Volume  fraction 


Figure  9.  Dependence  of  kt  on  the  volume 
fraction  of  the  single  crystal  in  the  1-3 
composite 


affect  the  input  mechanical  energy  density  in  the  energy  harvesting  cycle  using 


external  electric  boundary  condition  to  maximize  the  energy  output. 


(iv)  Experimental  investigation  on  materials  and  electric  boundary  conditions  to 
achieve  high  energy  density  harvested  with  high  efficiency.  As  a  first  step,  we 

compared  the  energy  density 
harvested  using  1-3  composites 
with  piezoceramic  vs.  1-3 
composites  with  piezo-single 
crystal.  Shown  in  figure  10  is 
the  schematic  of  the  experiment 
set-up  where  the  value  of  the 
load  resistor  Rl  is  determined  by 
the  electric  impedance  matching 
condition  to  the  electroactive 
polymer  to  reach  the  highest 

Figure  10.  Schematic  of  the  test  set-up  where  the  energy  density  harvested  in  this 

mechanical  stress  is  provided  and  regulated  by  an  Instron  configuration.  It  should  be  noted 
Universal  Test  Machine 
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that  the  efficiency  and  energy  output  of  this  simple  circuit  are  lower  than  these  listed 
in  Table  I  (as  discussed  in  (i)). 


In  the  following,  we  derive  the  efficiency  and  load  condition  for  the  set-up  in  figure 
10.  Suppose  the  stress,  strain,  and  electrical  displacement  are  non-zero  only  along  the 
z-direction  (see  figure  7),  the  constitutive  relations  of  the  electroactive  material  are  in 
the  form  of 


(19) 


S  =  sT  +  dE 
D  =  dT  +  eE 

where  S  and  T  are  the  strain  and  stress  component,  and  D  and  E  are  the  electric 
displacement  and  electric  field,  respectively,  s,  d,  and  8  are  the  elastic  compliance, 
the  piezoelectric  constant,  and  the  dielectric  permittivity.  For  the  1-3  PMN- 
33PT/Epoxy  composite,  Eq  (19)  can  be  rewritten  as: 


l  =  sF  +  dV, 
Q  =  dF  +  CV 


(20) 


where  the  dimensional  change  1  is  related  to  the  strain  S  as  l=St,  the  force  F=TA,  the 

charge  Q=DA,  the  voltage  V=Et.  The  spring  constant  s  =  s—  ,  the  capacitance 

A 

A 

C  =  s—  and  d  =  d .  A  is  the  cross  section  area  and  t  is  thickness  of  the  sample. 


When  a  sinusoidal  mechanical  force  F  =  F0  sin(<yt)  is  applied  on  the  sample  and  a 

sinusoid  voltage  would  be  generated,  the  phase  of  which  is  to  be  determined: 
V  -VQ  sin  (cot  +  6) ,  then  the  charge  and  current  are: 

Q  =  dF0  s'm(ax)  +  CV0  sin(£Ot  +  6) 

dQ  -  <21) 

i  =  —  =  dF0Q)cos(ax )  +  CV0cocos(cot  +  6) 
at 

Finally,  after  some  calculation,  the  power  dissipation  of  (or  delivered  to)  the  resistor  is 

(22) 


IW-g- 


2  R  2(1  +  Q)2C2R2) 

When  the  dissipation  resistor  is  equal  with  the  impedance  of  the  sample  voltage  source, 
the  power  is  the  maximum  power  dissipated  in  the  resistor,  so  the  power  peak  is  at 

(23) 

At  such  a  resistive  load,  the  output  power  density  is 

_  (Qc12Fq 

Power peak  =  — ^  (24) 

Since  the  input  peak  mechanical  energy  is  2 F^s  and  the  converted  electrical  for  one 


2  772 


cycle  is 
k: 


nd  F( 
2  C 


.  The  ratio  of  two  can  be  related  to  the  electromechanical  coupling  factor 
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^?2F02/2C_^2 

2 F025  ’’  4C$ 

Shown  in  figure  11  are  the  experimental 
piezoceramic  PZT. 


Jtk1 
~~4~ 
data  for 


the  1-3  composites  with 


Figure  11.  The  output  electric  power  density  from  a  1-3  composite  with  piezo-ceramic  PZT  (a)  without  DC 
bias  field  (left)  and  (b)  with  a  5.5  kV/cm  DC  bias  field  (right).  The  upper  curves  are  for  the  increased  stress 
and  the  lower  curves  are  for  the  reduced  stress  during  the  cycle.  Apparently,  without  DC  bias  field,  the 
piezocomposite  depoled.  No  depoling  occurs  under  DC  bias  field  and  the  output  power  density  is  higher. 


Power  density  for  PMN-PT  1-3  DC=200V 

Har.  Power  density  using  1-3  single  crystal 

dimension=(8.1 8mm*4.3mm*0.43mm) 

98- 

(dimenslon=1cm*1cm*0.068cm)  in  88.9  MPa 

Stress  (MPa) 


Figure  12.  The  output  electric  power  density  from  thel-3  composite  with  ferroelectric  relaxor  PMN-PT 
single  crystals  (a)  as  a  function  of  the  applied  stress  signal  (left,  for  the  composite  with  37%  crystals)  and 
under  4.6  kV/cm  DC  bias  and  (b)  under  a  stress  of  88.9  MPa  as  a  function  of  the  load  resistor  (right). 
The  maximum  output  power  density  reaches  96  mW/cc  at  10  Mfi  load.  The  operation  frequency  is  4  Hz. 
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Figure  12  presents  the  experimental  data  for  the  1-3  composites  with  ferroelectric 
relaxor  single  crystal  PMN-PT.  Figure  12(a)  is  for  the  1-3  composite  with  37%  of 
single  crystal.  In  the  experiment,  a  DC  bias  of  4.6  kV/cm  was  applied  to  stabilize  the 
polarization.  Figure  12(b)  shows  how  the  output  power  density  varies  with  the  load 
resistor  under  a  stress  cycle  with  Tmax=89  MPa,  here  the  1-3  composite  contains  56%  of 
PMN-PT  single  crystals.  The  output  power  density  can  reach  96  mW/cc  (at  4  Hz) 
when  the  resistor  load  is  optimized.  Table  III  summarizes  these  experimental  results. 

As  shown  from  figures  11  and  12,  the  1-3  composites  with  piezo-single  crystal  can 
generate  two  times  output  power  density  and  show  an  efficiency  more  than  2  times 
higher  in  comparison  with  that  obtained  from  1-3  composites  with  piezo-PZT  ceramics. 
The  results  are  consistent  with  the  fact  that  the  efficiency  of  1-3  single  crystal 
composites  is  twice  of  that  of  1-3  piezoceramic  composites. 

No  further  experiments  were  conducted  related  to  the  power  output  with  other  power 
electronic  circuits,  for  example,  the  open  circuit  electric  boundary  condition,  to 
maximize  the  system  efficiency  as  discussed  in  (i). 
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IV.  Summary 

The  results  from  this  preliminary  investigation  show  that  there  is  a  great  potential  in 
utilizing  the  newly  developed  electroactive  materials  for  the  electric  energy  harvesting 
during  walking  as  well  as  from  other  mechanical  sources.  Furthermore,  by  working  with 
smart  electronics  to  tune  the  material  parameters  during  the  energy  harvesting  cycle,  one 
can  achieve  much  higher  energy  efficiency  in  comparison  with  the  electroactive  material 
itself. 
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